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ABSTRACT 

This study deals with the effects of slip and phase difference in a steady flow of an incompressible 
asymmetric rigid channel with permeable walls. It is assumed that the effect of fluid absorption 
through permeable walls is accounted by prescribing flux as a function of axial distance. The 
perturbation method is applied to linearize the non-linear governing equations by assuming the ratio 
of inlet width to wavelength to be small. Effects of the above parameters on the velocity profile, 
mean pressure drop and wall shear stress are studied in detail and explained graphically. 
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1 INTRODUCTION 

The study of the flow of viscous fluid in an asymmetric channel of varying cross section with 
permeable walls is much interested in recent years in view of its numerous applications in many 
physiological and engineering problems. Fluid flow in renal tubules was studied by many authors. 
Mathematical modeling of the flow in proximal renal tubule was first studied by Macey [16] where 
he considered the flow of an incompressible viscous fluid through a circular tube with a linear rate of 
reabsorption at the wall. Bulk flow in the proximal tubule decays exponentially with the axial 
distance was calculated by Kelman [5]. Then, Macey [17] used this condition to solve the equations of 
motion to find the average pressure drop. Marshall et.o/[7] and Palatt et al[ 12] studied the physical 
conditions existing at the rigid permeable wall instead of prescribing the flux /radial velocity at the 
wall. 

In all the above studies the researcher considered the renal tubule to be symmetry. But in general, 
renal tubules may not be symmetric throughout their length. A hydrodynamical aspect of an 
incompressible viscous fluid in a circular tube of varying cross-section with reabsorption at the wall is 
studied by Radhakrishnamacharya et o/[14]. Flow in rigid tubes of slowly varying cross-section with 
absorbing wall is studied by Peeyush Chandra and Krishna Prasad [13]. Fluid flow through a 
diverging/converging tube with variable wall permeability was studied by Chaturani and Ranganatha 
[ 2 ]. 
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The concept of slowly varying flow is given by Manton [6] where he obtained an asymptotic series 
solution for the low Reynolds number flow through an axisymmetric tube, where radius varies slowly 
in the axial direction. 

The effects of slope parameter and reabsorption coefficient on the flow of fluid in a symmetric 
channel with varying cross section with no-slip velocity at the walls are studied by Muthu and 
Tesfahun [9]. 

In all the above studies the researchers have taken the boundary condition at the wall to be a no-slip 
condition, whereas the no-slip condition is one of the aspects on which the mechanics of the viscous 
liquids is built. However, there are many situations where this assumption does not hold [15]. 
Elshahed [8] illustrated the significance of the effect of slip at the wall. Also, the slip would be most 
useful for certain problems in chemical engineering and other applications ([15], [3], [4], [18], [19]). 
Fluid flow through the non-uniform channel with permeable wall and slip effect in symmetry channel 
is studied by [11]. Further, Muthu and Tesfahun [10] discussed the flow through in renal tubule by 
considering the asymmetric channel of varying cross-section, whereas Waseem et al [20] study the 
effect of slip on fluid flow in a channel of the slowly varying cross section. 

Thus, in this paper, an attempt is made to understand the flow through renal tubule of asymmetry 
channel of varying cross-section and a slip velocity at the walls of the channel. 

2 MATHEMATICAL FORMULATION 


Here we consider an incompressible fluid flow through the asymmetric channel with a slowly varying 
cross-section. The boundaries of the channel wall are taken by Muthu et o/[10]as 
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Where d x and d 2 are the half width of the channel from the x-axis to 7/j(x) and r/ 2 (x ) 
respectively at the inlet (at x=0), a, and b x are amplitudes and A is the wavelength further a x ,b l , 
d x , d 2 ,(f) satisfies the condition 
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Figure 1. Geometry of the channel 
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We shall consider the motion of the fluid to be laminar and steady and the channel to be long 
enough to neglect the initial and end effects. The equations of continuity and momentum are given 
by 
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Where u and v are the velocity components along the x andjy axes respectively, p is the 

u 

pressure, p is the density of the fluid and v = — is kinematic viscosity. 

P 

In order to complete the formulation of the problem, the boundary conditions are taken as follows, 
(a) The tangential velocity at the wall is not zero. That is. 
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Where P is slip parameter and y is the specific permeability of the porous medium. 

(b) The reabsorption has been accounted for by considering the bulk flow as a decreasing 
function of x. That is, the flux across a cross-section is given by 




Q(x)= J u(x,y)dy = Q o F(ax ) , 


( 8 ) 


Where F(ax ) = 1 when a = 0 and decreases with x, a > 0 is the reabsorption coefficient and is a 
constant, and Q a is the flux across the cross-section at x=0. 


The boundary conditions (6) and (7) are well known Beavers and Joseph[l] condition when applied to 
tangential velocity. 


We introduce the stream function y such that 
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Where d =d l +d 2 . 


By introducing the above non-dimensional variables the equations (3)— (5) can be written as (the 
primes are dropped) 
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Further, the boundary conditions (6—8) becomes 
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The parameter R is the Reynolds number and 8 is the wavenumber (the ratio of inlet width to the 
wavelength), e x and£ 2 are amplitude ratios (the ratios of amplitudes a x and b x to the inlet width 

respectively) and and P 2 are ratios of distance from the x-axis to the upper wall and lower wall to 

the inlet width respectively. In this problem, we consider exponentially decaying bulk flow [6] that is, 
in equation (8), F is taken as 

F(ax) = e- ax (15) 

3 METHOD OF SOLUTION 

It is observed that the flow is quite complex because of nonlinearity of governing equation and the 
boundary conditions (10)-(14). Thus to solve equation (10) for velocity components, in the present 
analysis, we assume the wave number SU (long wavelength approximation). We shall seek a 
solution for stream function \|/(x,y) in the form of a power series in terms of 8 , as 

V|/(x,j)=\|/ 0 (x,j) + ^\)/ 1 (x,j) + .. (16) 


Substituting equation (16) in equations (10)-(14), and equating the coefficients of like powers of 8 , 
we get the following sets of equations for v|/ 0 (x,y), \|/j (x,y),... 
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Zeroth Order System 


^ = 0 
&y 4 

The corresponding boundary conditions are: 

^ = t ;; = /7 0) and y = rj 2 (x ) 

dy dy 

\)/ 0 = F(ax ) = e~ ax at y = q, (x) 

\|/ 0 = - F(ax ) = - ^ e~“* a? j = r| 2 (x) 


(17) 

(18) 

(19) 

( 20 ) 


First Order System 


gVi 

5/ 


= R 


ay 0 ay 0 ay 0 

<9j 2 ck dx dy 3 


( 21 ) 


The corresponding boundary conditions are: 
d\\i 


= v4sin(2^x) 

dy y ’ dx l dy 

= Bsm{ 27 tx + (j)) 

dy dx 

Vo =° aty = rii(x) 


^-A^lnx) 8 ^ 


dydx 


^-B S in(2xx + t) S -^ 


dy 1 

and 7 = ri 2 (x) 


dydx 


at y = 7, ( x ) 

at y = rj(x) 


( 22 ) 

(23) 

(24) 


Similar expressions can be written for higher orders of 8 . However, since we are looking for an 
approximate analytical solution for the problem, we consider up to the order of S l equations. 


The solution of equation (17) along with the corresponding boundary conditions (18-20) as 
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Following the similar procedure as in equation (25) the solution of equation (21) along with boundary 
conditions (22-24) is 
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By substituting the value of v|/ 0 and v|/, in equation (16), we get 
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Now, the nondimensional pressure p(x,y^ can be obtained by using equations (27), (9) and (4), 
and it is given as 
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The mean pressure is given as 
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Further, the mean pressure drop between x = 0 and x = x 0 is 

&P( x o) = P(°)-P( x o) 

The wall shear stress r w (x) is defined as 
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Using the non-dimensional quantity T m = — r w andr W2 = , the wall shear stress and 
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It may be noted that in equation (28), the integrals are difficult to evaluate analytically to get closed 
form expression for p(x,y^j . Therefore, they are calculated by numerical integration. 


4 RESULTS AND DISCUSSION 

The purpose of the present discussion is to analyze the behavior of a steady incompressible fluid flow 
in an asymmetric channel of slowly varying cross-section with absorbing wall by considering a slip 
velocity at the walls. 
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It may be noted that (f) characterized the phase difference which varies in the range 0 < (f> < n . Here 
<j) = 0 represents symmetric channel. ^ = 7i represents the asymmetric channel with waves are in 
phase. Ctand^ represents reabsorption coefficient and slip at the channel wall, respectively. It is 
observed that in the absence of slip i.e., § — »■ 0 , our results are in tune with those of Muthu and 
Tesfahun[10]. 

We discuss the effect of these parameters on the transverse velocity v(x,y ) , mean pressure drop 
(A p) and wall shear stress (j w ). The following parameters are fixed as A=-0.0628, B=0.0628, 

f 

(3j = 0.5, p 2 = — 0.5, 8j = 0.1, s 2 = —0.1, 8 = 0.1 in our numerical calculation. For low Reynolds 
number flow, we have takeni? = 1.0. To see the effect of £ we have taken ^=0,0.15 and 0.4. 

The Transverse velocity v: 

The transverse velocity v(x,y ) which is obtained from equations (9) and (27) Here we have 
discussed the effects of the phase difference ( (j ) ), in the presence of non-zero slip coefficient (^ ) on 
the transverse velocity by taking the behavior at a different cross-section of the channel. We have 

taken x = 0.1, 0.5, 0.9 and^ = 0,^-,zr . 

Figure 2(a) displays the effect of ( (f>) on v at x=0.1 and ^=0.0. It may be observed that as ( (/>) 
increases from 0 to n , the magnitude of v decreases. It may be remarked that the reabsorption value 
at the wall is fixed at x=0.1 and when ( (j > ) increases, the cross-sectional area is reduced. This results 
in lesser v values. Now, if ^=0.15 similar effect is observed as above. When ( <f>) varies from 0 to n 
mixed trends is observed in velocity. 

If \ =0.4 the velocity decreases when (f> varies from 0 to K . But comparing with no slip ( \ =0.0) case 
velocity increases in quantity. This may be due to the effect of the slip (see figures 2(b), 2(c)). 

Figure 3(a) displays the effect of ^ on v at x=0.5 with ^=0.0. It may be noted that as 0 increases 
the magnitude of v has mixed trends, due to the variation of the cross-section of the channel at x=0.5. 

If ^=0.15 and ^=0.4, similar mixed trends is observed on v, due to the effect of slip (see figures 3(b), 
3(c)). Figures 4(a)-4(c) display the effect of ^ on v when x=0.9 for \ =0.0, 0.15 and 0.4. It is observed 
that as (j) increases the magnitude of v has a mixed trend. 

Mean Pressure drop 

The value of the mean pressure drop (29) over the length of the channel is calculated from different 
values of ^and^. Figure 5(a) represents the effect of (j> when ^=0.0. It is observed that as the 
width of channel contracts, the mean pressure drop increases. Particularly, at the entrance of the 
channel, the mean pressure drop for the asymmetrical channel is more than the symmetrical channel. 
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It can be understood from figure 5(a) for (f> = 0, — , n . However, due to contraction in the middle 
of the channel, the reverse is true at the end of the channel. 

When ^=0.15, a similar trend as mentioned above is observed, with a quality difference (see fig.5 
(b)). 

As t, =0.4, the trend is reversed, this shows the effect of slip combined with asymmetry nature of the 
channel (fig. 5( c) ). 

Magnitude of wall shear stress 

The effects of ^ and i;on the magnitude of the wall shear stress (| r | and| t u , 2 \ ) are presented in 
figures 6 and 7 respectively. 

It may be noted from figures 6(a) to 6(b), and 7(a) to 7(b) that the upper wall and lower on shear 
stress (in magnitude) increases as the channel changed from symmetry to asymmetry and no-slip to 
slip conditions, except in the middle of the channel where there is more contraction. 

This indicates that as the width of the channel decreases due to asymmetry nature of walls, 
(| r | and | r w | ) increases. But when ^=0.4, the nature of the curve is oscillatory [See fig. 6(c) and 7 


The main contribution of this study is to see the effect of the phase difference in the presence of slip 
at the walls on the flow of incompressible fluid in an asymmetric channel of the slowly varying cross- 
section. The mathematical problem is solved using a regular perturbation method assuming the ratio 
of inlet width to wavelength is small. We observe the following observation in the present study. 


(ii) As the channel changes from symmetric to asymmetric the mean pressure drop increases. 

(iii) The wall shear stress increases as the channel changes from symmetry to asymmetry and 
no slip to slip. 

APPENDIX - IMPORTANT FORMULAS AND GRAPHS 
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5 CONCLUSIONS 


(i) As phase difference increases the magnitude of velocity decreases. 
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Distribution of Transverse velocity v with y at 
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Figures 2(a)- 2(c) Distribution of Transverse velocity v with y atx=0.1, a=1.0 
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Distribution of Transverse velocity v with y at 
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Figures 3(a)- 3(c) Distribution of Transverse velocity v with y at x=0.5, a=1.0 
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Figures 4(a)- 4(c) Distribution of Transverse velocity v with y atx=0.9, a=1.0 
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Distribution of mean pressure drop (AP)with x. 



1 


9 = ° 9- - JT ■ ■ ■ 

■ ■ 6 = jt 


Distribution of mean pressure drop (AP)with x. 
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Distribution of mean pressure drop (SP)with x. 
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Figures 5(a)- 5(c) Distribution of Mean Pressure Drop with x 
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Distribution of magnitude of |Twl| with x =0 
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Distribution of magnitude of |Twl| withx 
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Distribution of magnitude of |Twl| with x 
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Figures 6(a)- 6(c) Distribution of Magnitude / Twlj of with x 
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Distribution of magnitude of |Tw2| with x 
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Distribution of magnitude of |Tw2| with x 
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Distribution of magnitude of |Tw2| with x 
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Figures 7(a)-7(c) Distribution of Magnitude / Tw2j of with x 
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